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ABSTRACT 

Current models of galaxy formation predict that galaxy pairs of comparable mag- 
nitudes should become increasingly rare with decreasing luminosity. This seems at 
odds with the relatively high frequency of pairings among dwarf galaxies in the Local 
Group. We use literature data to show that ~ 30% of all satellites of the Milky Way 
and Andromeda galaxies brighter than My = — 8 are found in likely physical pairs of 
comparable luminosity. Besides the previously recognised pairings of the Magellanic 
Clouds and of NGC 147/NGC 185, other candidate pairs include the Ursa Minor and 
Draco dwarf spheroidals, as well as the And I/And III satellites of M31. These pairs are 
much closer than expected by chance if the radial and angular distributions of satel- 
lites were uncorrelated; in addition, they have very similar line-of-sight velocities and 
luminosities that differ by less than three magnitudes. In contrast, the same criteria 
pair fewer than 4% of satellites in N-body/semi- analytic models that match the radial 
distribution and luminosity function of Local Group satellites. If confirmed in studies 
of larger samples, the high frequency of dwarf galaxy pairings may provide interesting 
clues to the formation of faint galaxies in the current cosmological paradigm. 

Key words: galaxies: dwarf - galaxies: formation - galaxies: evolution - Local Group 



1 INTRODUCTION 

The shapes of the galaxy and dark halo mass functions 
differ substantially in the ACDM paradigm (see, e.g., 
I Benson et alj [2003 ). This is usually interpreted to imply 
that the "efficiency" of galaxy formation, as measured by 
the ratio between the stellar mass of a galaxy (M ga i) and 
the viria{3 mass of its host halo (M200), varies strongly with 
virial mass. In particular, M gal /M2oo should decrease steeply 
toward low halo masses in order to match the shallow faint 
end of the galaxy luminosity function (see, e .g., lMoster et al.l 
l20ld ; lBehroozi et al J 1201(3 ; iGuo et al.|[201lb . 

On the scale of dwarf galaxies, which we define con- 
ventionally here as those with M ga i < 1O 9,5 M0, simple 
abundance-matching models suggest a dependenc e nearly as 
steep as M ga i tx M| 00 in the dwarf galaxy regime (|Guo et al.l 
l2010l h Such steep scaling would imply that dwarfs span- 
ning several decades in stellar mass should nevertheless in- 
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1 Wc define all virial quantities as those corresponding to a sphere 
of mean density equal to 200 times the critical density for clo- 
sure. Furthermore, we shall hereafter refer to the stellar mass of 
a galaxy as "galaxy mass" , for brevity. 



habit halos of similar virial mass. In addition, extrapolating 
such models to the faintest galaxies known indicate that 
few, if any, galaxies more massive than a few million solar 
masses are expected to form in halos with virial mass below 
10 10 M o . 

Recent work has highlighted potential disagreements 
between these model predictions and observations, includ- 
ing the lack of a c haracteristic veloci ty at the faint-end of 
blind HI surveys (|Zwaan et al.l l201ut ) ; and the low virial 
mass (substantially below 10 10 Mq) inferred from dynamical 
data for the dwarf spheroida l companions of the Milky Way 
(iBovlan-Kolchin et alj|2012f ) and for nearby dwarf irregulars 
^Ferrero et alj|2012j h It could be argued, however, that the 
evidence for substantial disagreement is unconvincing, given 
that the inferences are either indirect (in the case of the HI 
velocity function) or based on small and h eterogeneous sam- 
ples (in the case of th e nearby dwarfs; see lWang et aljfeoi^ : 
IVera-Ciro et~alll2012h . 

It is therefore important to consider further tests of 
the model predictions. We explore here how the steep M ga i- 
M200 relation predicted for dwarfs affects the frequency of 
galaxy pairs of comparable luminosity. Such pairs, when 
close enough to inhabit the same dark matter halo (referred 
to hereafter as "physical pairs"), are expected to be rare at 
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all luminosities, but especially so in the scale of dwarfs. This 
is because the fainter companions in physical pairs trace the 
halo substructure, and subhalos are, by and large, far less 
massive than the main halo: the most massive subhalo typi- 
cally has a mass only one hundredth that of the m ain system 
(see, e.g.. lSpringel et al.ll20o4 IWang et al.ll2012T l. 

Pairs of comparable luminosity are therefore more likely 
to form in fairly massive halos, where galaxy formation ef- 
ficiency decreases with increasing halo mass, partly com- 
pensating the mass difference between the main halo and 
its most massive subhalo. On dwarf galaxy scales the situ- 
ation is reversed, and the precipitous decline in galaxy for- 
mation efficiency with decreasing halo mass should curb the 
formation of physical pairs of comparable luminosity. More 
generally speaking, isolated associations of dw arf galaxies 
should be rare. Th ey are known to exist (e.g., iTullv et alj 
2006; Soa resl [20071 ) . but their cosmological abundance and 
dependence o n luminosity have not yet been adequately es- 
tablished fsee lSales et alj|2012l . for a recent attempt). 

The Local Group offers an interesting environment to 
test these ideas. Advantages include the fact that, away from 
the "zone of avoidance" caused by Galactic dust, the cen- 
sus of Milky Way (M W) satellites brighte r than M v ~ -8 
is complete (see, e.g., I Whiting et al.l [20071 ). and that accu- 
rate magnitudes, positions, distances, and line-of-sight ve- 
locities are known for all. Many fainter systems in the Lo- 
cal Group still remain undiscovered, as demonstrated by 
recent discoveries both by the Sloan D igital Sky Survey 
in the Milky Way (|Koposov et all 12008! ') . and by the Pan- 
Andromeda Archaeologic al Survey around M31 (PAndAS; 
iMcConnachie et al.llioog '). In the latter the census is likely 
complete within the PAndAS survey area down to a mag- 
nitude limit of Mv ~ —6.5, although some brighter dwarfs 
at larger radius likely have yet to be identified (as exem- 
plified by the d i scoveries of Andro meda XXVIII and XXIX; 
iBell et alJ[20Tll : ISlater et al.ll201ll ). 

Further, we know that at least some of the satellites of 
the Milky Way and M31 are very likely physically associ- 
ated and bound to each other. An obvious pairing is that 
of the Magellanic Clouds (see, e.g.. iKallivavalil et al.|[2006l . 
and references therein). Around M31, there have been sug- 
gestions that NGC 14 7 and NGC 185 also form a bound pair 
( van den Berghlll998l ). If these pairs are bound, their vicin- 
ity to their primary galaxy suggests that we are observing 
them just before they are separated by the tidal f ield of 
the main galaxy (|Besla et al.ll2007l ; ISales et alj|201lh . This 
implies very recent accretion and indicates that their occur- 
rence should not be uncommon amongst isolated systems. 

Taken at face value, the existence of these two pairs of 
dwarfs seems at odds with the expected rarity of such associ- 
ations. We use this as motivation to search, using literature 
data, for other dwarf galaxy pairs in the Local Group. We de- 
scribe in Sec.[2]the observational dataset and the simulated 
satellite dataset we use for comparison. In Sec. [3] we intro- 
duce the pairing procedure we have adopted and compare 
the results with those obtained when the same procedure 
is applied to a hybrid N-body/semi-analytic model of satel- 
lite galaxy formation applied to N-body simulations from 
the Aquarius Project. We conclude with a brief summary in 
Seed 
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We use the recent compilation by IMcConnachie] i|2012f ) as 
the source of the positions, distances, line-of-sight veloci- 
ties, and magnitudes of Local Group dwarfs that we use in 
our analysis. All velocities are heliocentric and corrected to 
the rest frame of the Galaxy. In order to prevent biases due 
to incompleteness, we consider only satellites brighter than 
M v = -8 located within 300 kpc of the Milky Way or An- 
dromeda (M31) galaxies. The sample consists of 29 dwarfs, 
17 of which orbit around M31; the rest are satellites of the 
Milky Way. 

For comparison, we have identified analogous samples 
of simulated satellites in the six ~ 10 12 M Q halos of the 
Aquarius Project dSpringel et all 12008!) using the model of 
IStarkenburg et all 121 '). This is a semi-analytic model 
grafted onto the level-2 Aquarius runs, which simulate each 
halo with several hundred million particles, thus ensuring a 
resolution high enough to track the formation of all halos 
and subhalos that might plausibly host the dwarf galaxies 
brighter than Mv = — 8 we use in our analysis. The model 
satellites of each Aquarius halo have luminosity and radial 
distributions that are broadly consistent with the Milky Way 
and M31 and therefore provide a useful testbed of the sig- 
nificance of our results for ACDM dwarf galaxy formation 
models. 

There are a total of 175 simulated satellites brighter 
than Mv = —8 within 300 kpc of the primary galaxies of 
all six Aquarius halos (on average 29 per halo). The model 
provides not only the full 3D position and velocity infor- 
mation for all of them, but also allows us to track their 
evolution. Our simulated sample does not include satellites 
whose dark matter halos have been fully disrupted by tides, 
si nce their fate is un c ertain . We refer the interested reader 
to lStarkenburg et al.l (|2012l ) for details on the semi-analytic 
model. 



3 ANALYSIS AND RESULTS 

It has long been noticed that the spatial di stribution of 
Milky Way satellites is highly anisotropic (|Lvnden-Belll 
Il976h . and is often described as a polar plane whose signifi- 
cance has been the matter of much recent debate (see, e.g., 
Kroupa et all 120051; [Zentner et all 120051 ; iLibeskind et all 



20051 ; iMetz et al.ll2007h . Around M31. 14 out of 17 satel- 



lites i n our sample are in the hem isphere nearer the Milky 
Way l|McConnachie &i Irw in 2006). Further, several have re- 
cently been shown to delineate a flattened structure that 
in total comprises at least half of all known M31 satellites 
ijlbata et al.l 120131 ). These are unlikely configurations for a 
virialized population and hint strongly at recent accretion. 

Our pairing procedure begins by identifying satellites 
whose nearest neighbour is unusually close when compared 
with the probability distribution of nearest-neighbour dis- 
tances, d nn , obtained by Monte Carlo sampling a random 
isotropic population of satellites with the same total number 
and radial distribution. We illustrate this in Fig.[T] where we 
show the d nn distribution expected for two satellites of the 
Milky Way (Draco and Sagittarius), and two of M31 (And 
III and And IX). 

The bottom left panel of Fig. [T] shows that the near- 
est satellite to Sagittarius (the LMC, 52 kpc away) is about 



Dra And II 




Dist. to Near. Neighb. [kpc] Dist. to Near. Neighb. [kpc] 

Figure 1. Distribution of nearest-neighbour distances, d n n, to 
two Milky Way satellites (panels on the left) and two M31 satel- 
lites (panels on the right) expected if satellites were distributed 
isotropically about each primary with a radial distribution consis- 
tent with the observed one. The probability. P, that a satellite's 
nearest neighbour lies, by chance, as close as or closer than ob- 
served is highlighted by the shaded region of each histogram and 
quoted in each panel's legend. A downward arrow indicates the 
distance to the primary galaxy. The top panels illustrate cases 
where the probability is rather small, indicative of a potential 
physical association. The bottom panels, on the other hand, illus- 
trate two cases where the nearest neighbours are not significantly 
closer than expected at random. 



twice as far as the distance at which the probability distribu- 
tion of nearest-neighbour distances peaks. If the radial and 
angular distribution of Milky Way satellites were uncorre- 
lated then Sagittarius would be expected to have a nearest 
neighbour as close or closer than observed in 99 out of 100 
random realizations (P — 0.99). Sagittarius is thus relatively 
isolated and unlikely to be a member of a physical pair. 

The situation reverses for Draco: its nearest neighbour, 
Ursa Minor, lies only 23 kpc away. This is much closer than 
expected at random; a nearest neighbour that close occurs in 
fewer than 4 out of 100 random trials (P = 0.04). The right- 
hand panels of Fig. [1] show as well two analogous examples 
for the M31 satellite population. In this case, And IX is 
unlikely to be a member of a pair (P = 0.90), whereas And 
III is unusually close to And I (P = 0.10), hinting at a 
possible physical association. 

The distribution of the probability, P, that d nn is as 
small or smaller than observed if the radial and angular dis- 
tribution of Milky Way satellites were uncorrelated is shown 
by the solid and dotted lines in Fig. [2] for the Milky Way 
and M31 satellites, respectively. Both curves are rather sim- 
ilar, indicating that the P distribution is insensitive to the 
total number of satellites or their radial distribution. It is 
also insensitive to assuming that the satellite distribution 
is isotropic. Indeed, the solid and dotted curves in Fig. [2] 
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log (P) 

Figure 2. Distribution of the probability, P, of having a nearest 
neighbour as close as or closer than observed if the satellites were 
isotropically distributed around each primary and had the same 
radial distribution as that of the Milky Way (solid thick line) 
and M31 (solid dashed line). Only 20% of satellites are expected 
to have P < 0.2, with a very weak dependence on the number 
of satellites and the shape of the radial profile. The probability 
distribution obtained for a semi-analytic model applied to the 
six Galaxy-sized halos of the Aquarius Project is shown by the 
dashed blue histogram; 48 out of 175 satellites have P < 0.2, or 
27% of the total. The corresponding distribution for Local Group 
satellites is shown by the solid red histogram (the contribution 
of Milky Way satellites is highlighted by the shaded area of the 
histogram). In the Local Group, more than 40% of satellites have 
P < 0.2, a result expected to happen by chance in fewer than one 
in 100 random realizations. The angular and radial distributions 
of satellites thus seem highly correlated and suggest the presence 
of physically-associated pairs. 



change almost imperceptibly if we confine the Monte Carlo 
samples to a three-dimensional structure as flat as observed 
for the Milky Way, i.e., roughly 3:1 in its major-to-minor 
axis ratio. 

On the other hand, these distributions differ markedly 
from the results for the Milky Way (shaded histogram in 
Fig. [H or the combined M31+MW satellites (labelled "Lo- 
cal Group" in Fig. [2j . A K-S test yields a probability of less 
than 0.3% that the Local Group P distribution is statisti- 
cally consistent with that of the random samples. There is a 
clear excess of smaller-than-expected nearest-neighbour dis- 
tances in the Local Group that is difficult to account just 
by chance. For example, 45% of Local Group satellites have 
P < 0.2 compared with the 20% expected if the distribution 
was isotropic. 

Our pairing procedure therefore retains all P < 0.2 
pairs (listed horizontally in the labels of Fig. [2]| for further 
scrutiny. A true physical pair must also differ little in ve- 
locity, so we impose a maximum difference of 75 km/s in 
the line-of-sight velocity difference of the likely members. 
This threshold is motivated by the velocity difference of the 
Magellanic Clouds, where there is little doubt about their 
physical association. We assume for simplicity that the same 
threshold applies regardless of the luminosity of the pair; 
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Figure 3. Left panel: Galactocontric velocity versus distance for Local Group satellites brighter than My = —8. Distances are measured 
from the center of each primary: MW satellites are shown as magenta squares; M31's as green triangles. Velocities for the former 
are Galactocontric radial velocities; for the latter they refer to line-of-sight velocities relative to the systemic velocity of M31. Dotted 
curves indicate, for reference, the escape velocity from an NFW halo with virial velocity V200 = 250 km/s and concentration c = 10 
l|Navarro et al.lll997f ). Filled symbols highlight satellites with a nearest neighbour much closer than expected by chance (P < 0.2, see 
Fig. [2}- Pairs satisfying additional proximity criteria in velocity (AV < 75 km/s) and magnitude (AMy < 3) are joined together by 
ellipses to indicate that they arc likely physical pairs. These constitute 28% of the total and include (i) the Magellanic Clouds; (ii) NGC 
147 and NGC 185; (iii) Ursa Minor and Draco; and (iv) And I and And III. Right panel: Same as left panel, but for the semi-analytic 
satellite population of the six Aquarius halos. Different colors correspond to different halos. Note that the same criteria that pair 28% of 
Local Group satellites link only six satellites in the Aquarius simulations, or just 3% of the total. 



this is also consistent with the idea that most d warfs should 
inhab it halos of similar mass (see Sec. [TJ (see I Sales et al.l 
12011 for a more thorough discussion). Finally, since we are 
mainly interested in pairs of comparable luminosities we re- 
tain as likely pairs only those where their magnitudes differ 
by AM V < 3. 

Four pairs remain after applying these constraints, as 
shown in Fig. [3] Around the Milky Way, aside from the 
Magellanic Clouds, the Ursa Minor and Draco pair is singled 
out: they are 23 kpc apart, and their velocities differ by 
only 12 km/s. Further, they have both approximately the 
same luminosity, which makes them especially singular. The 
procedure also joins together two pairs of satellites around 
M31: NGC 147 and NGC185, as well as And I and And III. 
The latter are separated by 33 kpc; differ in velocity by just 
32 km/s; and in magnitude by 1.7. 

This analysis suggests that nearly 30% of Local Group 
satellites are in likely pairs (8 out of 29). Of the four pairs, 
two are almost indisputably associated (the Magellanic 
Clou ds and NGC 147/NGC 185; see, however, iGeha et al.l 
|2010| for an alternate view of the latter) but the other two 
might in principle result from chance close encounters be- 
tween unrelated satellites where projection effects reduce the 
line-of-sight velocity difference. In order to quantify these 
effects we have applied the same pairing procedure to the 
satellite populations of the si x Aquarius halos, a s ident ified 
by the semi-analytic model of IStarkenburg et al.l |20ll . 



The blue dotted histogram in Fig. [2] shows that the P 
distribution for Aquarius satellites differs little from that 
expected from an isotropic distribution. Only 27% of the 
175 satellites have P < 0.2; of those only 3 pairs (i.e., fewer 
than 4%) pass as well the velocity and magnitude criteria. 
The three Aquarius pairs (out of six halos) singled out by 
the analysis are shown in the right-hand panel of Fig. [3] 
Tracking their orbits back in time reveals that none of them 
are actually physically related but that they result simply 
from chance, transient associations in position and veloc- 
ity space. Note that this does not imply that all satellites 
have b een accreted in isolation. As discussed bv lWang et al.l 
|20ll . a few of the bright satellites in Aquarius were ac- 
creted in groups, but the accretion happened early and the 
groups have long been disrupted by the tidal field of the 
main halo. This confirms the model expectation that physi- 
cal associations amongst satellites should be extremely rare. 
The relatively high frequency of satellite pairings in the Lo- 
cal Group indicates that the radial and angular distributions 
of satellites are correlated, a fact that is not easily accounted 
for by current dwarf galaxy formation models in the ACDM 
paradigm. 
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4 SUMMARY AND CONCLUSIONS 

We have studied possible pairings amongst the satellites of 
the Milky Way and of M31. Our procedure, which iden- 
tifies unusually close associations in position and velocity 
space, suggests that 8 out of the 29 satellites brighter than 
Mv = —8 (i.e., nearly ~ 30%) form 4 likely pairs of compa- 
rable luminosity (AMy < 3). These include the Magellanic 
Clouds; Ursa Minor and Draco; NGC 147 and NGC 185; as 
well as And I and And III. 

The same pairing procedure applied to a semi-analytic 
model of the satellite population in the six halos of the 
Aquarius Project yields a likely pair fraction of fewer than 
4%, even though the model satellites have luminosity and ra- 
dial distributions that match closely that of the Local Group 
spirals. As expected, none of the Aquarius pairs correspond 
to true binary systems; rather, they result from transient as- 
sociations between otherwise unrelated satellites. The high 
pair frequency of the Local Group is unlikely to be just a 
statistical fluke: the likely pair fraction of Aquarius satellites 
never exceeds 12% in thousands of random trials where their 
magnitudes, angular directions and velocities are reshuffled. 

We interpret these results as indicative of signifi- 
cant clustering in the dwarf galaxy population of the Lo- 
cal Group. Although our analysis only considers satel- 
lites brighter than Mv = —8 due to incompleteness con- 
cerns and in order to allow comparison with simulations, 
there have also been suggestions that some of the fainter 
Galactic satellit e s are found in associations. In particular, 
iBelokurov et alj (|2008T ) show that Leo IV and Leo V are 
close to each other spatially and differ little in their line- 
of-sight velocities. Further, the association of galaxies clos- 
est to the Local Gr oup, sometimes referr ed to as the loose 
NGC 3109 group (|van den Berghl 1 19991 ). consists of four 
dwarf galaxies (NGC3109, Antlia, Sextans A and B), two of 
which appear to be int eracting (NGC3019 and Antlia; see 
iBarnes fc de Blokll200lf ). These associations are expected to 
be rare according to current ACDM galaxy formation mod- 
els, and may signal the presence of a mechanism that boosts 
the likelihood of forming dwarf galaxies near other dwarfs, 
favouring some environments over others. No such effect is 
present in current models of dwarf galaxy formation, which 
rely largely on the mass accretion history of individual halos 
to set the properties of a dwarf. 

Candidate mechanisms that may bias the regions where 
dwarf galaxies form in a way that enh ances their clusterin g 
include patchy reionization (see, e.g.. iLunnan et al]|2012h ; 
large-scal e feedback effects on the surroundings of massive 
galax ies (| Alvarez et al.ll2009l ; iBusha et al.ll2010l ; iFont et all 
l201ll ) ; and the interaction of dwarf galaxies with the cosmic 
web (|Bem'tez-Llambav et al.ll2012l ). However, the effects of 
these models on dwarf galaxies have yet to be developed 
fully and as a consequence the importance of such effects on 
the dwarf galaxy population at large is still unknown. 

It is therefore important to firm up these findings (i) 
by extending the analysis to fainter satellites, which should 
be possible once photometric surveys of the northern and 
southern sky extend the complete catalog of Milky Way 
satellites to fainter magnitudes; (ii) by verifying, through 
accurate proper motion studies, that the associations in po- 
sition and line-of-sight velocities remain once the full 6D 
phase space information is considered; (iii) by searching 



for relic evidence of past interactions between likely pairs 
(such as the Magellanic Stream for the LMC/ SMC (see, e.g., 
iMathewson et all 1 1974 iPutman et"all ll99ST ); and, finally, 
(iv) by extending this kind of analysis to a volume-limited 
surv ey of dwarf galaxy associations in the local universe (see, 
e.g., iKarachentsev fc Makarovll2008l ). If confirmed, the en- 
hanced clustering of dwarfs may offer important clues to the 
formation of faint galaxies that have yet to be identified and 
fully incorporated into galaxy formation models. 
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